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ABSTRACT 



Context. It has recently been proposed that the jets of low-luminosity radio galaxies are powered by direct accretion of the hot phase 
of the IGM onto the central black hole. Cold gas remains a plausible alternative fuel supply, however. The most compelling evidence 
that cold gas plays a role in fueling radio galaxies is that dust is detected more commonly and/or in larger quantities in (elliptical) 
radio galaxies compared with radio-quiet elliptical galaxies. On the other hand, only small numbers of radio galaxies have yet been 
detected in CO (and even fewer imaged), and whether or not all radio galaxies have enough cold gas to fuel their jets remains an 
open question. If so, then the dynamics of the cold gas in the nuclei of radio galaxies may provide important clues to the fuelling 
mechanism. 

Aims. The only instrument capable of imaging the molecular component on scales relevant to the accretion process is ALMA, but 
very little is yet known about CO in southern radio galaxies. Our aim is to measure the CO content in a complete volume-limited 
sample of southern radio galaxies, in order to create a well-defined list of nearby targets to be imaged in the near future with ALMA. 
Methods. APEX** has recently been equipped with a receiver (APEX-1) able to observe the 230 GHz waveband. This allows us to 
search for CO(2 - 1) line emission in our target galaxies. 

Results. Here we present the results for our first three southern targets, proposed for APEX- 1 spectroscopy during science verification: 
NGC 3557, IC4296 and NGC 1399. The experiment was successful with two targets detected, and possible indications for a double- 
homed CO line profile, consistent with ordered rotation. These early results are encouraging, demonstrating that APEX can efficiently 
detect CO in nearby radio galaxies. We therefore plan to use APEX to obtain CO spectroscopy for all our southern targets. 
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1. Introduction 



It i s gene rally thought (see e.g. Heckman et al. ( 1986); Baum et 
al. 0992)) that powerful radio galaxies (typically Fanaroff-Riley 
type II or FR II, Fanaroff & Riley ( 19741) ) are triggered by merg- 
ers or collisions with gas-rich galaxies. During the merger event, 
part of the gas can be transported to the central region of the 
merging system, where it may trigger a starburst, and, if it loses 
sufficient angular momentum, feed the AGN. The radio jets may 
then be powered either from the gravitational energy released 
during the accretion process, or by electromagnetic extract ion of 
the spin energy of the black hole (Blandford & Znajek (1 19771) 1. 
which also requires an accretion flow to anchor a magnetic field. 

On the other hand, it has been suggested recently that accre- 
tion in low luminosity (or more generally low excitation) radio 
galaxies, often, but not exclusively, Fanaroff-Riley type I or FR I, 
may occur dir ectly f rom the hot phas e of th e interstellar medium 
(Allen et al. (120061) : Evans et al. (120081) ). Hot (T > 10 7 K) 
plasma is observed to be ubiquitous in radio-galaxy nuclei, and 
Evans et al. ((2008)) suggest that the gravitational energy re- 
leased if accretion occurs at the Bondi rate may be sufficient to 
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** This publication is based on data acquired with the Atacama 
Pathfinder Experiment (APEX). APEX is a collaboration between 
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Observatory, and the Onsala Space Observatory. 



power the jets in FRI radio galaxies. It is unclear whether this 
rate is relevant, however, since it assumes spherically symmetric 
accretion and takes no account of the need for angular momen- 
tum loss. The accretion rate in the vicinity of the black hole can- 
not be estimated directly and may be much lower. In any case, 
McNamara et al. ((2010)) find that accretion at the Bondi rate 
is generally unable to power the jets in brightest cluster galaxies 
and argue that electromagnetic extraction of the black hole spin 
energy is required. 



An alternative fuel source is cold gas, which is often present 
in large amounts, and is detectable via emission and absorption 
by dust, HI and CO. This cold gas quite often comes in the form 
of dusty disks, as ca n be c learly seen in high resolution HST 
images (Capetti et al. (120001) : Verdoes-Kleijn (1 19991) : de Ruiter et 
al. (|2002)). Dust is observed in 53% of the B2 sample of nearby 
radio galaxies (mostly FR I) a nd the dust mass is correlated with 
radio power (de Ruiter et al. (2002)). Similar results have been 
found for a sample of nearby 3CR radio galaxies. At least 30% 
of the sources show evidence for dust absorption. Dust in FR I 
radio galaxies is generally situated in well defined disks on small 
(< 2.5 kpc) scales, and the radio source axis tends to lie nearly 
perpendicular to the dust disk. Very interestingly FR I galaxies 
have derived dust masses that are typically larger than the dust 
masses fou nd in a matched sample of radio-quiet ellipticals (de 
Koff et al. (120001) ). If the mass of dust and cold gas are also 
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correlated, as seems plausible, then this association argues that 
accretion of cold gas may indeed power the radio jets. 

If accretion of (initially) cold gas is the dominant fuelling 
mechanism, then we should be able to detect the gas in every 
low-excitation radio galaxy. In normal, extended, FR I sources 
very little HI gas is found close to the nucleus (certainly no 
more than in early-type galaxies of similar mass selected with- 
out reference to radio emission); in most cases only upper limits 
M ffl ^ 10 8 M o are available (Emonts et al. (120101) '). On the other 
hand, we know from several observing programmes that CO ob- 
servations of nearby radio galaxies produce a high detection rate 
(see e.g. Pr andon i et al. d2007l) : Ocana Flaquer et al. d2010b '). 
Lim et al. ( d2000l) ) detected 12 CO(l - 0) and 12 CO(2 - 1) emis- 
sion from the FR I radio galaxies 3C 3 1 and 264 with the IRAM 
30m Telescope and established that the line profiles indicate 
disk r otatio n; this appears to be a common property (Prandoni 
et al. d2007h 1. 

The derived masses of molecular gas are certainly adequate 
to power the jets in at least some cases. For example, the well- 
studied FRI radio galaxy 3C31 h as an estimated jet power 
Pjet ~ lO^ergs 1 (Laing & Bridle d2002l) ). while the lifetime 
of the source must be > 10 8 /(v/0.01c) years, where v is the ad- 
vance velocity of the tails. This suggests that the lifetime of 3C3 1 
may be appreciably longer than those of FRII sources (10 6 - 10 7 
years; Blundell & Rawlings (2000)) and therefore that a higher 
molecular gas mass may be needed to fuel it, despite its lower 
radio luminosity. For an assumed conversion efficiency rj from 
the rate of gravitational energy release to jet power, the requir oH 
mass accretion rate for 3C 31 is M — Pjet/^c 2 w O.O2M yr~ 
r] = 0. 1 and the molecular gas mass needed to fuel the AGN 
10 8 years would be ~ 2 x 10 6 MO (assuming a constant acc 
tion rate). The mass of mo lecular gas associated with 3C31 t 
~ 10 9 M o (Lim et al. d2000l) V This implies that, even if the e j 
ciency is significantly smaller or the advance speed is < 0.0 
there is ample fuel for any reasonable source lifetime. 3C 
is a luminous FRI radio galaxy, and correspondingly smal 
amounts of cold gas would be required to fuel less power 
sources, assuming efficient accretion. Current Hi mass upp^ 
limits (of the order of 5 X 10 7 - 10 8 MO; e.g. Fig. [3]) do not rule 
out the possibility that cool gas is the dominant fuelling mecha- 
nism in low luminosity radio galaxies, as well as in their brighter 
counterparts. If this is the case, we expect to detect CO in virtu- 
ally all radio galaxies, provided that CO spectroscopy goes deep 
enough. 

Only high-resolution imaging of CO or other molecular 
lines would allow us to probe the dynamics of the cool gas di- 
rectly, to determine if and how much of it is infalling, and what 
kind of motions (rotation, non-circular, radial in- or out-flows) 
are present. Interferometry at mm wavelengths is challenging 
with current instrumentation, and has been possible for only a 
few objec ts. Interferometric observations of 3C 31 by Okuda et 
al. ((12005)) showed that the CO coin cides spatially with the dust 
disk observed by HST (Martel et al. d 19991) ) and is in ordered ro- 
tation. These authors suggest that the cold gas is in stable orbits, 
in which case the actual accretion rate may be very low. 

ALMA, thanks to its combination of angular resolution and 
sensitivity, will give us a unique opportunity to assess the cold 
accretion scenario in low luminosity radio galaxies, by directly 
imaging the gas on scales ~10 - lOOpc and determining its dy- 
namical state. Very little is known about cold gas in Southern 
galaxies, however, with very few objects observed in CO. 

In this paper we report the results of pilot APEX- 1 CO line 
observations of Southern radio galaxies. In Sect. 2 the Southern 
radio galaxy sample is presented with particular emphasis on the 




Fig. 1. 1.4 GHz VLA radio continuum contours (left) and HST 
images (right) for radio galaxies NGC 1399 (top), NGC3557 
(middle), IC4296 (bottom). 1.4 GHz radio continuum contour 
images are taken respectively from Shurkin et al. ((2008)) for 
NGC 1399 (here overlaid to a Chand ra X-r ay color image); from 
the NVSS database ( Condo n et al. dl998l) ) for NGC 3557; and 
from Killeen et al. (( fl986h) fo r IC4296. The HST images are 
taken from Lauer et al. ((2005|)). 

three sources discussed in this paper. In Sect. 3 the APEX line 
measurements are described, while in Sect. 4 a brief discussion 
of the results is reported. In the following we assume a standard 
A-CDM cosmology with H = 70 kms 1 Mpc" 1 . 

2. Sample Description 

The primary radio-galaxy samples (3C, B2) for which CO obser- 
vations have been made are difficult to observe from the ALMA 
site. We have therefore built a volume limited sample of low lu- 
minosity radio galaxies in the Southern sky, for which we intend 
to obtain complete single-dish CO coverage in order to prepare 
for ALMA imaging. This sample has been defined in a very sim- 
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Table 1. APEX-1 CO(2 - 1) line measurements 



Source 
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(min) 
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PKS0336-35 


[NGC 1399] 


0.0048 


108.2 


272 


40 
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0.20 


127 


2.1 


7.47 
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0.35 


0.39 


365 


3.3 


7.75 


PKS 1107-372 


[NGC 3557] 


0.0102 


18.4 


316 


40 


2 


1.52 


248 


3.5 


9.02 


PKS1333-33 


[IC4296] 


0.0125 


37.4 


315 


40 


1 


< 0.51 






< 8.70 



ilar way to the Northern sample we extracted from the B2 radio- 
source catalogue and foll owed up in CO with the IRAM 30m 
telescope (Prandoni et al. d2007l) ). 

The parent sample is selected from the Parkes 2.7-GHz sur- 
vey, in the decl ination range -17° < 6 < -40°, as described 
by Ekers et al. ( (fl98l ). It has a radio flux-density limit of 0.25 
Jy at 2.7 GHz and an optical limit of my < 17.0, and consists 
of 191 radio galaxies. We extracted those sources with redshifts 
Z < 0.03 which are associated with elliptical galaxies. The re- 
sulting 11 objects thus form a small, but volume-limited and 
complete subset of the Ekers et al. (( 1989)) sample. All eleven 
sources of the volume-limited sample are of the FRI type. The 
subset is quite similar to the sample of 3C and B2 radio galaxies 
in the Northern Hemisphere, which consists of 23 low luminos- 
ity radio sources associated with elliptical galaxies at z < 0.03; 
16 of these have been observed in CO with the IRAM 30m tele- 
scope (Prandoni et al. (2007)). 

As part of the science verification programme of the APEX- 
1 receiver, three of the eleven sources in our subsample 
were selected for observation in the 12 CO(2 - 1) line. They 
were: PKS0336-35 (NGC 1399), PKS1107-372 (NGC3557), 
and PKS 1 333-33 (IC 4296). All three are radio sources with twin 
radio jets (see Fig. [TJ left panels), and had been obser ved wi th 
the WPFC2 of the Hubble Space Telescope (Lauer et al. (120051) ). 
While no trace of dust was found in the nucleus of NGC 1399, 
the other two show prominent dusty disks around their respective 
nuclei (see Fig.Q] right panels). 

A fourth sample member, PKS0320-37 (Fornax A) was im- 
aged in 12 CO(l -0) and 12 CO(2- 1) by Horellou et al. ((2003)); it 
shows a complex dust structure in HST observations (Grillmair 
et al. (fl986k 

No CO imaging is available for the remaining sample mem- 
bers. 



3. Line Observations and Measurements 

We used the APEX-1 SHeFI receiver connected to the FFTS 
backend to search for emission in the (2-1) transition of 12 CO. 
Low luminosity radio galaxies have 12 CO(2- 1) line width s typ i- 
cally spanning 200 -600 kms 1 FWHM (Prandoni et al. ((2007)) 
and are therefore fully sampled by a 1-GHz wide individual unit 
of the FFTS backend (covering a total velocity range of 1300 
kms 1 at 230 GHz). The 30 arcsec FWHM beam probes the 
presence of molecular gas in the galaxy cores (the inner 3-7 
kpc, depending on the redshift), allowing a direct comparison 
with dust structures on similar scales imaged with HST. 

The observations were carried out in Summer 2008 and the 
data were reduced with the CLASS package. Scans severely af- 
fected by ripples and/or noisy baselines were removed and first- 
order baselines were subtracted. The noise levels obtained were 
typically r™ s ^1-2 mK (lcr) at a reference smoothing width 
of Av * 40 kms -1 . Line fluxes were measured by numerically 



integrating over the channels in the line profile. Line widths were 
measured as full widths at half power. A source was considered 
detected when the 12 CO(2 - 1) emission line had T v a eak > 3r™\ 

and tentatively detected when T^ eak > 2T™ S . In case of non 
d etect ions, upper limits were calculated, following Evans et al. 
( (120051) ). through the relation 



T a Av (Kkms -1 ) < 



3T rms < Ay > 

V< Av > /Av re . 



(1) 



where < Av > is the mean FWHM line width for radio galax- 
ies (in ou r case assume d equ al to 500 kms -1 , see Prandoni et 
al. (120071) . Evans et al. (2005)); < Av res > is the velocity resolu- 
tion (in our case m 40 km s -1 after smoothing) and 



T a 



(2) 



where Av is the channel width (« 30 MHz, after smoothing) in 
Hz and t is the on-source integration time in seconds. 

H2 molecular masses were deriv ed using the same rela- 
tion^ used by Gordon et al. ((1992)) and later by Lim et al. 
( (120001) 1. but modified for a A-CDM cosmology and H = 70 
kms -1 Mpc -1 . We assumed an aperture efficiency rj a = 0.63 for 
230-GHz APEX observations. 

A summary of our CO line measurements is given in Table[TJ 
where, for each source, we list the source redshift (z); the on- 
source integration time (few), after flagging data affected by poor 
baselines and/or ripples; the mean system temperature during the 
observations (T sys ); the channel width after smoothing (Av res ) 
and the corresponding noise level in the line spectra (T™ s ); the 
line integrated flux ("£T a dv); the FWHM line width (Av FWH m); 
the line signal to noise (Tf ak IT r a ms ) and the molecular gas mass 
(logM//,). The final CO line spectra are shown in Fig. [2] 

NGC 3557 was easily detected at a channel resolution of 
Av res ~ 40 km s _1 , while for IC4296 (the source at highest 
redshift) only an upper limit could be derived. In the case of 
NGC 1399 only the higher of the two peaks in the spectrum 
(Fig. |21 top left) has a signal-to-noise ratio T% e /T™ s > 2 and 
can be considered tentatively detected. However the systemic 
velocity of the galaxy falls between the two peak. By increas- 
ing the smoothing to Av res ~ 120 km s -1 (Fig. |2l top right), 
we can reduce the noise level to T' a ' ns ~ 0.35 mK and get a 
clear > 3<x detection over a wider velocity range. In this case 
we obtain I,T a dv = 0.39 K km s _1 , Av FWlffl = 365 kms -1 , and 
log M H2 1 MO = 7.75 (see TableQ}- NGC 1399 does not show de- 
tectable dust features in the HST image, and the presence of CO 
is particularly interesting in this case. There is some evidence 
for a double-horned structure in the spectrum of NGC 3557, al- 
though the statistical significance is low (Fig. [2] bottom left). If 
confirmed, the presence of CO in ordered rotation would be con- 
sistent with the strong dusty disk seen in HST observations (see 
Fig.[TJ. A very marked dusty disk is present even in IC4296 (see 
Fig.[TJ and we suspect that a longer observation (toN 30 min), 
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Fig. 2. 12 CO(2 - 1) spectra obtained at APEX for the radio galaxies NGC 1399 (top), NGC 3557 (bottom, left) and IC 4296 (bottom, 
right). NGC 1399 is shown for two different smoothing boxes: a velocity width of » 40 kms _I (top, left), as for the other two 
sources; and a wider velocity width of * 120 kms~' (top, right). Typical spectrum error bars (±T™ 1S and ±Av res /2) are shown for 
reference at the systemic velocity of the galaxies. Horizontal dashed lines indicate the zero level T a = 0. 



possibly with better weather conditions (T sys < 300 K), will lead 
to a detection. 



4. Discussion and Conclusions 

The H2 masses of the four low-luminosity southern radio galax- 
ies observed so far (including Fornax A) are plotted against red- 
shift as magenta symbols in Fig. [3j together with the 23 B 2 
z < 0.03 radio galaxies (blue symbols, Prandoni et al. (120071) ') . 
t he bri ghter z < 0.03 3C radio galaxies studied by Lim et al. 
( (l2003h : green symbols) and the z < 0.0 233 UGC galaxies with 
radio jets studied by Leon et al. ((2003); red symbols). We note 
that the three northern samples are partially overlapping and that 
some upper limits have have been updated using recent mea- 
surements by Ocana Flacquer et al. ((2010)). The deepest CO 
line observations with existing facilities have 3<x detection limits 
=s 3 x 0.1 mK, as indicated by the solid lines in Fig. [3] CO de- 
tection rates vary between 45% and 65% for the three samples, 
but these differences are not significant once account is taken of 
the relative sensitivities and redshift distributions. We find that 
molecular gas masses are very similar over the full range of red- 
shift (and, implicitly, radio power) for the combined samples, 
spanning the range * 10 7 — 10 9 M , with upper limits varying 



from « 10 7 M to w 10 8 M Q , depending on distance. These 
values are consistent with the hypothesis that the jets of low- 
luminosity (FRI) radio galaxies can be powered by accretion 
of cold gas, but there is a large scatter in molecular gas mass 
a t a giv en radio luminosity, as also found by McNamara et al. 
( (120101) 1 for brightest cluster galaxies. 

Seven examples of CO detections without visible dust, in- 
cluding NGC 1399, are plotted in Fig. [3] The inferred H2 masse s 
range from 8 x 10 7 to 4 x 10 s Mo (Ocana Flaquer et al. (2010); 
this paper). It is therefore not yet clear whether there is a sig- 
nificant difference between the H2 masses of radio galaxies with 
and without dust. Nevertheless CO line detections are frequent 
in radio galaxies with dust, even more when the dust is in form 
of nuclear disks. In the sample of 23 e lliptical B2 radio galaxies 
with z < 0.03 (Prandoni et al. (|2007|)), seventeen objects have 
detailed optical imaging (mostly with HST), and ten show clear 
dusty disks. CO line observations are available for nine of these 
objects, and m 80% (7/9) were detected. We also find that al- 
most one half of the objects with dusty disks detected in CO have 
double-horned line profiles. The double-horned CO line profiles 
characteristic of ordered rotation (indicated by the crosses in 
Figj3j occur over a wide range of molecular gas masses. Since 
many of the other detections are still marginal, it is possible that 
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Fig. 3. H2 mass vs redshift for the z < 0.03 B2 radio galax- 
ies (large blue empty symbols); for z < 0.03 3C radio galax- 
ies (small green empty symbols); for z < 0.0233 UGC galaxies 
(small red filled symbols); and for the four (including Fornax A) 
z < 0.03 Southern radio galaxies observed so far (large magenta 
filled symbols). Circles refer to detections and triangles to upper 
limits. Crosses indicate sources which display a double-horned 
CO line profile. Solid lines indicate the 3cr detection limit 
for c (= T™ s ) = 0.1, 0.2 and 0.4 mK respectively (assuming a 
500 kms" 1 wide CO line, see equation[TJ. The vertical dashed 
line at z — 0.0233 indicates the redshift limit to which all three 
Northern samples overlap. 



the molecular gas in the majority of low-luminosity radio galax- 
ies is in the form of kpc-scale nuclear disks or rings, cospatial 
with the dust and in ordered rotation. There are exceptions, how- 
ever: Fornax A represents a more complex case where the CO 
shows a multi-component line profile associated with the com- 
plex dust structures visible in HST B - 1 images. The most prob- 
able scenario is that the gas comes from the ac cretio n of one or 
more small gas-rich galaxies (Horellou et al. J200P). The fact 
that we are finding CO in radio galaxies with no detectable dust 
may suggest that in (at least some) objects the bulk of the gas is 
in a cold (molecular) phase, supporting the idea that virtually all 
radio galaxies may be fuelled by cold gas. 

These results raise a number of questions which can be an- 
swered by observations at optical, mid-infrared and mm wave- 
lengths, and in particular by sensitive, high-resolution imaging 
with ALMA, as follows. 

1. Is the molecular gas in disks or in rings with well-defined 
inner edges? 

2. Is the majority of the gas in stable orbits? 

3. What is the evidence for non-circular motions, especially in- 
fall? 

4. Where does the gas come from: is its angular momentum 
consistent with that of the stellar population or different, in- 
dicating an external origin? 



5. Do the molecular and ionized gas distributions have the same 
spatial distributions and kinematics? Does the gas become 
progressively hotter and then ionized closer to the nucleus? 

6. Is there evidence for ongoing or recently-completed star for- 
mation? 

In conclusion, we have shown that observations with APEX 
in the 230-GHz band can efficiently detect CO in low-luminosity 
radio galaxies. We plan to use APEX to obtain CO spectroscopy 
for the entire southern sample and to build up a target list for 
imaging with ALMA. 

Acknowledgements. The authors thank the APEX team, which successfully 
overcame a number of technical obstacles to obtain the APEX- 1 science verifi- 
cation observations presented here. The authors would like to thank the referee, 
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